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(57) Abstract: A method for the adaptive adjustment of a PMD compensator in optical fiber communication systems with the com- 
pensator comprising a cascade of adjustable optical devices through which passes an optical signal to be compensated and comprising 
the steps of computing the Stokes parameters So, S lt S 2 , S 3 in a number Q of different frequencies of the signal output from the com- 
pensator, producing control signals for parameters of at least some of said adjustable optical devices so as to make virtually constant 
said Stokes parameters computed at different frequencies. A compensator comprising a cascade of adjustable optical devices (12-14) 
through which passes an optical signal to be compensated, an adjustment system which takes the components yi(t) e y 2 (t) on the two 
orthogonal polarizations from the signal at the compensator output, and which comprises a controller (15, 16) which on the basis of 
said components computes the Stokes parameters S 0 , S u S 2 , S 3 in a number Q of different frequencies of the signal output by the 
compensator and which emits control signals for at least some of said adjustable optical devices so as to make virtually constant the 
Stokes parameters computed at the different frequencies. 
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POLARISATION MODE DISPERSION COMPENSATOR 



The present invention relates to methods of adaptive adjustment of PMD compensators 
in optical fiber communication systems. The present invention also relates to a 
compensator in accordance with said method. 

5 In optical fiber telecommunications equipment the need to compensate the effects of 
polarization mode dispersion (PMD) which occur when an optical signal travels in an 
optical fiber based connection is known. 

It is known that PMD causes distortion and dispersion of optical signals sent over 
10 optical fiber connections making the signals distorted and dispersed. The different time 
delays among the various signal components in the various polarization states acquire 
increasing importance with the increase in transmission speeds. In modern optical fiber 
based transmission systems with ever higher frequencies (10 Gbit/s and more), accurate 
compensation of PMD effects becomes very important and delicate. This compensation 
15 must be dynamic and performed at adequate speed. 

The general purpose of the present invention is to remedy the above mentioned 
shortcomings by making available a method of fast, accurate adaptive adjustment of a 
PMD compensator and a compensator in accordance with said method. 
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In view of this purpose it was sought to provide in accordance with the present 
invention a method for the adaptive adjustment of a PMD compensator in optical fiber 
communication systems with the compensator comprising a cascade of adjustable 
optical devices over which passes an optical signal to be compensated comprising the 
5 steps of computing the Stokes parameters SO, SI, S2, S3 in a number Q of different 
frequencies of the signal output from the compensator, producing control signals for 
parameters of at least some of said adjustable optical devices so as to make virtually 
constant said Stokes parameters computed at the different frequencies. 

10 In accordance with the present invention it was also sought to realize a PMD 
compensator in optical fiber communication systems applying the method and 
comprising a cascade of adjustable optical devices over which passes an optical signal 
to be compensated and an adjustment system which takes the components yi(t) and y 2 (t) 
on the two orthogonal polarizations at the compensator output with the adjustment 

15 system comprising a controller which on the basis of said components taken computes 
the Stokes parameters So, Si, S 2 , S 3 in a number Q of different frequencies of the signal 
output from the compensator and which emits control signals for at least some of said 
adjustable optical devices so as to make virtually constant the Stokes parameters 
computed at the different frequencies. 

20 

To clarify the explanation of the innovative principles of the present invention and its 
advantages compared with the prior art there is described below with the aid of the 
annexed drawings a possible embodiment thereof by way of non-limiting example 
applying said principles. In the drawings - 
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FIG 1 shows a block diagram of a PMD compensator with associated control circuit, 
and 

FIG 2 shows an equivalent model of the PMD compensator. 

5 

With reference to the FIGS FIG 1 shows the structure of a PMD compensator 
designated as a whole by reference number 10. This structure consists of the cascade of 
some optical devices which receive the signal from the transmission fiber 11. The first 
optical device is a polarization controller 12 (PC) which allows modification of the 
10 optical signal polarization at its input. There are three polarization maintaining fibers 
13 (PMF) separated by two optical rotators 14. 

A PMF fiber is a fiber which introduces a predetermined differential unit delay (DGD) 
between the components of the optical signal on the two principal states of polarization 
15 (PSP) termed slow PSP and fast PSP. 

In the case of the compensator shown in FIG 1 the DGD delays at the frequency of the 
optical carrier introduced by the three PMFs are respectively t c , oct c and (1-cc) t c with 
0<cc<l and with x c and a which are design parameters. 



20 



An optical rotator is a device which can change the polarization of the optical signal 
upon its input by an angle Q t (the figure shows for the first rotator and 9 2 for the 
second) on a maximum circle on the Poincarfc sphere. 
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An optical rotator is implemented in practice by means of a properly controlled PC. 

In FIG 1, xi(t) and x 2 (t) designate the components on the two PSPs of the optical signal 
at the compensator input whereas similarly yi(t) and y 2 (t) are the components of the 
5 optical signal at the compensator output. 



The input-output behavior of each optical device is described here by means of the so 
called Jones transfer matrix H(co) which is a 2 x 2 matrix characterized by frequency 
dependent components. Designating by W t (cfl) e W 2 (co) the Fourier transforms of the 
10 optical signal components at the device input the Fourier transforms Zi(co) e Z2(co) of 
the optical signal components at the device output are given by: 



(z { m 

Z 2 {co) 



= H«y) 



(1) 



15 



Thus the Jones transfer matrix of the PC is: 



(2) 



where hi e h 2 satisfy the condition |hi| 2 +|h 2 | 2 =1 and are frequency independent. 
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Denoting by <|>i and the PC control angles, hi and h 2 are expressed by: 
h!=-cos(4»2-<t)i)+ j sin(<|>2-<t>i) sin«J>i 
h 2 =-j sinOfe-^i) cos<|»i 



(3) 
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Clearly if the PC is controlled using other angles or voltages, different relationships will 
correlate these other parameters with hi and h2. The straightforward changes in the 
algorithms for adaptive adjustment of the PMD compensator are discussed below. 



Similarly, an optical rotator with rotation angle 9* is characterized by the following 
Jones matrix: 

f cos 0, sin 0, ^ (4) 
-sin0, cos0, 



10 The Jones transfer matrix of a PMF with DGD X\ may be expressed as RDR" 1 where D 
is defined as: 



D 



( e J<*r t n 



0 e 



0 



(5) 



15 and R is a unitary rotation matrix accounting for the PSPs' orientation. This matrix R 
may be taken as the identity matrix I without loss of generality when the PSPs of all the 
PMFs are aligned. 



20 



As shown in FIG 1, to control the PMD compensator a controller 15 is needed to 
produce optical device control signals of the compensator computed on the basis of the 
quantities sent to it by a controller pilot 16 termed controller driver (CD). 
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The CD feeds the controller with the quantities needed to update the compensator 
optical device control parameters. As described below, these quantities will be 
extracted by the CD from the signals at the input and/or output of the compensator. 

5 The controller will operate following the criterion described below and will use one of 
the two algorithms described below. 

To illustrate the PMD compensator adaptive adjustment algorithms let us assume that 
the controller can directly control the parameters fc v fc, 0i and 9 2 which we consolidate 
10 in a vector 6 defined as: 



If it is not so, in general there will be other parameters to control, for example some 
voltages, which will be linked to the previous ones in known relationships. 

15 

The time instants in which the update of the compensator parameters is realized are 
designated t n (con n=0,l,2...,), and T u designates the time interval between two 
successive updates, thus tn + i=t n +T u . In addition, Q(Q designates the value of the 
compensator parameters after the nth update. 

20 

In accordance with the method of the present invention the criterion for adjusting the 
compensator parameters employs the so-called Stokes parameters. Computation of the 
Stokes parameters for an optical signal is well known to those skilled in the art and is 
not further described. 
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Again in accordance with the method the parameters 9 of the compensator are adjusted 
to make constant the Stokes parameters computed at different frequencies on the 
compensator output signal. The four Stokes parameters So, S u S 2 e S 3 computed at the 
frequency ft are designated by: 

S Q \ fafl - S 0 j 
s i I/-/, - s u 

S 2 I/-/, " S 2J 

10 Similarly, the Stokes parameters computed at the frequency f p are designated by S 0(P , 
Si, p , S2, P eS3, p . 



Using these Stokes parameters the following unitary vectors are constructed with 
components given by the three Stokes parameters Si, S 2 , S 3 normalized at the parameter 
15 S 0 . (.) T below designates the transpose while (./ designates the complex conjugate: 



u 



and 



^2J ^xi 



'o,/ 
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\p 



>£p_ 



In the absence of PMD these two vectors are parallel. Consequently, if their quadratic 
Euclidean distance is considered Gip(8): 
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3 U 



>0,P 



* 5 V 

^2,1 + 



5,, 



J 0,p 



(6) 



5 which is a function of the parameters 0 of the PMD compensator it will be zero when 
the PMD is compensated at the two frequencies considered f/ and f p . 



Now consider a number Q of frequencies f,, /=1,2,...,Q. Compute the Stokes parameters 
at these frequencies and construct the corresponding units defined as explained above, 
10 i.e. with components given by the three Stokes parameters Si, S2, S3 normalized with 
respect to the parameter S 0 . All these units are parallel if and only if the sum of their 
quadratic Euclidean distances is zero. 

Consequently, to adaptively adjust the PMD compensator parameters we define the 
15 function G(0) which is to be minimized as the sum of the quadratic distances G !p (0) 
with l,p =1,2,...,Q, i.e. the sum of the quadratic distances of the pair of vectors at the 
different frequencies f/ and f p , for J,p=l,2,...Q: 

G(6) = ix<V 0 ) (7) 

20 

The update rule for the compensator parameters to be used in accordance with the 
present invention are: 
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3G(9) 



e = 6(r„) 



3G(9) 



^2 



e=e(U 



/=2 p=l »Pl 
/=2 p=l #2 



e=eoj 



e = 6(r n ) 



(8) 



dG(Q) 



3G(0) 



e=e(o 



30, 



e = e(o 



/«2 p=l 9^2 



10 where y>0 is a scale factor which controls the amount of the adjustment. 

In vector notation this means that the vector of the compensator parameters is updated 
by adding a new vector with its norm proportionate to the norm of the gradient of G(0) 
and with opposite direction, i.e. with all its components having their sign changed. This 
15 way, we are sure to move towards a relative minimum of the function G(9). 

All this is equivalent to: 



Q l-l 



(9) 
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A simplified version of (9) consists of an update by means of a constant norm vector 
and therefore an update which uses only the information on the direction of VG(9) . In 
this case the update rule becomes. 



i \ 
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e<U) = «('.) Vign VG(9)l 




= 6(0 ysign 



1=2 pn\ 



(10) 



10 



15 



where sign (z) designates a vector with unitary components and of the same sign as the 
components or the vector z. 

Two methods are now described for computing the gradient of the G(0) function and 
obtaining the required control parameters. 

First Method 

To implement the update rules (8) the partial derivatives of G(0) for 0 = 9 (to) can be 
computed using the following five-step procedure. 

Step 1. find the value of G[9(tn)]=G[<t)i(t n ), (^(tn), 0i(tn), 9 2 (t n )] at iteration n. To 
do this, in the time interval (t n , tn+Tu/5) the Stokes parameters at the above mentioned Q 
frequencies are derived and the value of the function G(9) is computed using equations 
(6) and (7) . 

Step 2, find the partial derivative 



3G(0) 



9 <*i 0 = e(f„) 



at iteration n. To do this, parameter fy\ is set at <t>i(t n )+A while the other 
parameters are left unchanged. The corresponding value of G(9), i.e. 
G[<t>i(t n )+A, <l>2(tn), 9i(t n ), 9 2 (t n )L is computed as in step 1 but in the time interval 
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(t n +T tt /5, t n +2Tu/5). The estimate of the partial derivative of G(9) as a function of 
<J>i is computed as: 



3G(9) 



e=e<o 



_ Gfe(» l ,) + A,# > (>.),g l (t.),g a (r.)]-Gk(».).^(r.),g 1 (ag > fr.)1 (11) 

A 



Step 3. Find the partial derivative: 



3G(6) 



9 = 9(0 



10 



15 



at iteration n. To do this the parameter is set at <fe(tn)+A while the other 
parameters are left changed. The corresponding value of G(8), i.e. G[<j>i(t„), 
MO+A, 9i(t n ), 62(tn)], )], is computed as in step 1 but in the time interval 
(t„+2Tu/5, t„+3T,/5). The estimate of the partial derivative of G(9) with respect 
to fa is computed as: 



3G(8) 



6 = 9(rj 



A 

(12) 
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Step 4: Find the partial derivative: 

3G(9) 



dd i 9 = 8(r B ) 

at iteration n. To do this, parameter 9i is set at 9i(tn)+ A while the other 
parameters are left unchanged, the corresponding value of G(9), i.e. G[<j>i(tn), 
<t»2(to), 9i(tn)+A, 9 2 (t n )], is computed as in Step 1 but in the time interval 
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(t n +3T u /5, t n +4Tu/5) and the estimate of the partial derivative of G(8) with 
respect to 0i is computed as: 



0G(0) 



00, 



e=0(O 



C?k(f a )>0 2 a B ).g i qj + A > <? 2 (fj]-Gk(r i ,),^a M ),g 1 a B ),g 2 a n )] 

A (13) 



Step 5: Find the partial derivative: 
0G(0) 



d0 2 e = e(o 

at iteration n. To do this the parameter & is set at MW+A while the other parameters 
10 are left changed. The corresponding value of G(0), i.e. GE^t,,), (^(to), 02(t„) +A], 
is computed as in step 1 but in the time interval (t„+4T u /5, t n +T u ). The estimate of the 
partial derivative of G(0) with respect to is computed as: 



0G(0) 



00, 



Gka B ).^aj^i(^).^^)-«-A]-G^(f n ),^(f n ),g,a n ),0 2 (fj] 



2 9 = 6(0 



(14) 
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The above parameter update is done only after estimation of the gradient has been 
completed. 



Note that in this case it is not necessary that the relationship between the control 
20 parameters of PC and optical rotators and the corresponding Jones matrices be known. 



Indeed, the partial derivatives of the function with respect to the compensator control 
parameters are computed without knowledge of this relationship. Consequently if the 
control parameters are different from those assumed as an example and are for example 
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some voltage or some other angle, we may similarly compute the partial derivative and 
update these different control parameters accordingly. 



Lastly, it is noted that when this algorithm is used the CD must receive only the optical 
5 signal at the compensator output and must supply the controller with the Stokes 
parameters computed at the Q frequencies f/, /=1,2,..,,Q. 

Second method 

10 When an accurate characterization of the PC and of each optical rotator is available the 
update rules can be expressed as a function of the signals on the two orthogonal 
polarizations at the compensator input and output. 

In this case, for the sake of convenience it is best to avoid normalization of the three 
15 Stokes parameters Si, S 2 e S 3 with respect to S 0 and use the function H(6) defined as: 

/=2 p=l (15) 

where 

H/p(e)KS l ,rS 1 /4.(S 2|I .S 2t p) 2 +(S3,rS 3 ,p) 2 (16) 



20 



Consequently we have new update rules similar to those expressed by equation (8) or 
equivalently (9) with the only change being that the new function H(6) must substitute 
the previous G(0). 
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Before describing how the gradient of this new function H(0) is to be computed let us 
introduce for convenient an equivalent model of the PMD compensator. 

Indeed it was found that the PMD compensator shown in FIG 1 is equivalent to a two- 
5 dimensional transversal filter with four tapped delay lines (TDL) combining the signals 
on the two principal polarization states (PSP). This equivalent model is shown in FIG 2 
where: 



For the sake of convenience let c(0)designate the vector whose components are the ci in 
(17). It is noted that the tap coefficients q of the four TDLs are not independent of each 
other. On the contrary, given four of them the others are completely determined by 
(17). In the FIG for the sake of clarity it is designated P=l-a. 



10 



c x ^cos^cos^ 
c 2 =-smG l smd 2 h l 
c 3 ^-sin^cosdy^ 
c 4 = -cos0!sin0 2 /£ 
c s = cos 0j cos #2^2 
c 6 =- sin 0 { sin djv^ 
c 7 =sm6 l cos0 2 h^ 
c 8 ^cos^sin^/i,* 



(17) 
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The gradient of Hi p (0) with respect to 6 is to be computed as follows: 



VH „ (6) = 4(S U - S hp ) ReU- (/)a[ (r) - y u {t)b] (0 - y hp (t)a T p (r) + y %p (t)b T p (f )] dtj 
+ - 5 2iP ) Re|~f +1 \,l (0a[ (0 + y u »b[ « - y 2tP {tW p (0 - y[ p <f)b] (t)]dt3 
- 4(5 3 , - S 3tp ) Imj^f 1 fc, (0a[ (0 + y u (t)bj (t) - (t)a T p (r) - y^ (t)b T p (t)]dtf 
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where: 

y ; ,/(t) and y 2 ,K0 are the signals yi(t) and y 2 (t) at the compensator output 
respectively filtered through a narrow band filter centered on the frequency f/ (similarly 
for y; ;p (t) and y 2 , P (t)); 
5 - a/(t) and bj(t) are the vectors: 



10 



^(r-T c -arr c ) 

x 2 At-ar c ) 
^ 2J (t-T c -ar e ) J 



b,(0 = 



' < ( (r-2r c ) 

*« (»-*«> 

-x* u (t-2T c ) 
-x u (t-/3r c ) j 
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with x/,i(t) and x 2 ,/(0 which are respectively the signals xi(t) and x 2 (t) at the 
compensator input filtered by a narrow band filter centered on the frequency fj 
(similarly for y/, p (t) and y 2 , p (t)); 
- J is the Jacobean matrix of the transformation c=c(9) defined as 



20 



'dc, 










Wz 


de l 


w 2 


dc 2 


dc 2 


dc 2 


dc 2 


ajj" 


df 2 




dd 2 






3c 8 








de l 


de 2 

2 J 



(18) 



The parameters 9 are updated in accordance with the rule 



e (-i 



ee B+1 ) = ecf„)-r£2vj^<wL _ e . 



(19) 
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or in accordance with the following simplified rule based only on the sign: 

Q 1 1 . 

ee* + i)=e(0 ysign ^J^Mq^) (20) 

5 When the control parameters are different from those taken as examples we will 
naturally have different relationships between these control parameters and the 
coefficients Ci. 

For example, if the PC is controlled by means of some voltages, given the relationship 
10 between these voltages and the coefficients hi and h 2 which appear in (2), by using the 
equations (17) we will be able to express the coefficients Ci as a function of these new 
control parameters. 

Consequently in computing the gradient of the function H(6), the only change we have 
15 to allow for is the expression of the Jacobean matrix J, which has to be changed 
accordingly. 

Lastly it is noted that when this second method is used the CD must receive the optical 
signals at the input and output of the compensator. The CD must supply the controller 
20 not only with the Stokes parameters for the optical signal at the compensator output and 
computed at the Q frequencies f/, /=1,2,...,Q but also with the signals x; fi (t), x 2 ,/(t), y;.;(t) 
e y 2 /t) corresponding to the Q frequencies f/, /=1,2,...,Q. 
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It is now clear that the predetermined purposes have been achieved by making available 
an effective method for adaptive control of a PMD compensator and a compensator 
applying this method. 

5 Naturally the above description of an embodiment applying the innovative principles of 
the present invention is given by way of non-limiting example of said principles within 
the scope of the exclusive right claimed here. 
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CLAIMS 

1. Method for the adaptive adjustment of a PMD compensator in optical fiber 
communication systems with the compensator comprising a cascade of adjustable 
optical devices over which passes an optical signal to be compensated comprising the 
steps of: 

computing the Stokes parameters S 0 , S u S 2 , S 3 in a number Q of different 
frequencies of the compensator output signal, and 

producing control signals for parameters of at least some of said 
adjustable optical devices so as to make virtually constant said Stokes 
parameters computed at different frequencies. 

2. Method in accordance with claim 1 comprising the steps of computing the 
Stokes parameters in pairs of frequencies fl and fp with l,p=l,2,..»,Q, obtaining at the 1th 
and pth frequencies of the Q frequencies the two series of Stokes parameters So.i, S u , 
S2,i, S3,! and So, p , Si, p , S2, p , S3 tP , computing a vector function of each series of Stokes 
parameters and producing the control signals in such a manner that said vectors 
function of the two series of parameters are virtually parallel to each other. 

3. Method in accordance with claim 2 in which said vectors are unitary norm 
vectors with components given by the Stokes parameters Si, S 2 , S 3 normalized to the 
Stokes parameter S 0 , i.e.: 
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and 



J 0,p »Q,p 



Method in accordance with claim 3 in which is defined the function: 

Q m 



G(e)=£S G /p( 0 ) 

/«2 p=l 



with 



<ye>= 



5 0,/ S 0,p , 



and the control signals are produced to minimize said function G(0). 



5. Method in accordance with claim 4 in which the optical devices comprise a 
polarization controller with controllable angles fo, fa and two rotators with controllable 
rotation angles respectively 0t and 62, and to minimize the function G(0) the updating of 
<|>i, $2, 81 and 82 of the compensator observes the following rules to go from the nth 
iteration to the n+lth iteration: 



3G(9) 



0 = 0(0 l ° 2 pal 



(6) 



0 = 0(0 



dG(B) 



e=e(o 



Z=2 p=*l </^2 



e=0(O 



3G(6) 



0-6(0 



30, 



e=e(o 



/«2 p=l VV\ 



e = e(».) 



^2 p=l 9# 2 



o=e(o 
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6. Method in accordance with claim 5 in which the partial derivatives of G(0) for 

0=9(t n) are computed in accordance with the following steps: 

Step 1. find the value of G[0(t n )]=G[(|),(t I1 ), <t>2(t n ), 0,(t D ), 0 2 (t n )] at 

iteration n; to do this, in the time interval (t n , t n +T u /5) the Stokes parameters at 

the Q frequencies are derived and the value of the function G(0) is computed. 

Step 2. find the partial derivative 
dG{9) 



301 



0 = 6>(fJ 



at iteration n; to do this, parameter §\ is set at <|>i(t n )+A while the other 
parameters are left unchanged, the corresponding value of G(0), i.e. GtMO+A, 
teCn), 0i(tn), 02(t n )L is computed as in step 1 but in the time interval (t n +T u /5, 
t n +2T u /5) and the estimate of the partial derivative of G(0) with respect to §\ is 
computed as: 



3G(0) 



e = 0(O A 



Step 3. Find the partial derivative: 
dG(9) 



at iteration n; to do this the parameter (fe is set at MtnHA while the other 
parameters are left changed, the corresponding value of G(0), i.e. G[<>i(t n ), 
<t>2(t n )+A, 0i(t n ), 02(t n )], is computed as in step 1 but in the time interval 
(t n +2T u /5, t n +3T u /5) and the estimate of the partial derivative of G(0) with 
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respect to fa is computed as: 



8G(8) 



e = 6(o 



■ gfe(<.).^fr,)+A,g 1 (f.Xg a fr,)]-gfe(r,),# 1 (f,),g 1 (/.),g a (r,)] 

A 



Step 4. Find the partial derivative: 
3G(8) 



30, 



at iteration n; to do this, parameter 0! is set at 9i(t n )+ A while the other 
parameters are left unchanged, the corresponding value of G(8), i.e. G[^>i(t n ) 
•fcCtn), 6i(tn)+A, 9 2 (t n )], is computed as in Step 1 but in the time interval 
(ta+STJS, t n +4Tu/5) and the estimate of the partial derivative of G(0) with 
respect to 6i is computed as: 



8G(9) 



B6 l 



, G|> I (f n ),^(rJ,e i (fJ + A,g 2 (f n )3-G[i» 1 (^), ( » 2 (?J,g 1 (f a ),g 2 ( fB )] 



e = 6(f B ) 



- Step 5. Find the partial derivative: 

8G(9) 



de 2 



6 = 9(0 



at iteration n; to do this the parameter (fe is set at (J)2(tn)+A while the other 
parameters are left unchanged, the corresponding value of G(0), i.e. G[<|)i(t Q ), 
<fe(tn), Bi(t n ), 9 2 (t n ) +A], is computed as in step 1 but in the time interval 
(tn+4T u /5, tn+Tu) and the estimate of the partial derivative of G(0) with respect to 
fa is computed as: 



9G(9) 

Be, 



e=e(o 
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7. Method in accordance with claim 1 comprising the steps of computing the 
Stokes parameters in pairs of frequencies f/ and f p with /,p=l, 2,..,Q, to obtain at the 1th 
and pth frequencies of the Q frequencies the two series of Stokes parameters Si,/, S2,/, 
S3,/ e Sij, S2,p, S 3iP , defining the function: 

#(6) = £5X0) 

with ^p(0)KSi^i/+(S2^S2,p)VS3,i-S 3 ,p) 2 
and producing said control signals to minimize said function H(d). 

8. Method in accordance with claim 7 in which the optical devices comprise a 
polarization controller with controllable angles (j>i, fa and two rotators with controllable 
rotation angles respectively 6i and 8 2 , and for minimizing the function H(Q) the 

updating of <j>i, (j>2, 9i and 02 of the compensator follows the following rules for passing 
from the nth iteration to the n+lthe iteration: 



e(^) = e^)-r||v^(e)| 0=ean) 

or the following simplified rule: 



Q(t n+l ) = d(t n )-)si g n 



L tl™» (e) le=e<,„) 



with VH U (0) equal to the gradient of Hip(9) with respect to 9 
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9. Method in accordance with claim 1 in which the PMD compensator is modeled 
like a two-dimensional transversal filter with four tappered delay lines combining the 
signals on the two principal states of polarization (PSP). 

10. Method in accordance with claim 9 in which the gradient VH^ifi) with respect 
to 6 is computed as: 



where yy.Kt), Y2,/(t) and yy, p (t), y 2 , p (t) are respectively the components yi(t) e 
y2(t) on the two orthogonal polarizations of the compensator output signal 
filtered respectively through a narrow band filter centered on the frequency 
f/ and i p ; and 

a/(t) e b/(t) are the vectors: 




*u(0 
xj,,(r-ar c ) 

xu(f-r € -ar e ) 

xuit-ar c ) 
*2.,(r-T c ) 



\ 



xl(t-T c -0T c ) 



a,(f) = 



b,(0 = 



-<,(f-2r c ) 



J 
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with x /( /(t) and x 2 ,/(t) which are respectively signals x^t) and x 2 (t) on the two 
orthogonal polarizations of the compensator input signal filtered with a narrow 
band filter centered on the frequency f/ (similarly a/(t) and b/(t) for y JtP (t) and 
y2, P (t)) with the frequency f/), and 

J is the Jacobean matrix of the transformation c=c(0) defined as 



J = 



f dc x 


3q 


3cj 






W 2 




de 2 




dc 2 


dc 2 


dc 2 




df 2 


de t 


d? 2 












at 


d9, 


W 2J 



(18) 



with c !? . . ,.,c 8 which are the tap coefficients of the four tappered delay lines. 



11. Method in accordance with claim 7 in which said parameters are consolidated in 
a vector 0 which is updated in accordance with the rule 

ec»)=e(«.)-r||v H , ( ^ e = 9(0 (19) 

or the following simplified rule based only on the sign: 



e('„ +1 )=e(0-pfcn 



(20) 



with VH y, (0) equal to the gradient of Hi p (9) with respect to 0 
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12. Method in accordance with claim 1 in which said optical devices comprise a 
polarization controller with control angles fa, fa and two optical rotators with rotation 
angles 0i and 02 and said parameters comprise said control angles fa, fa and said 
rotation angles Gi, 02. 

13. Method in accordance with claim 11 in which between the controller and an 
optical rotator and between optical rotators there are fibers which introduce a 
predetermined differential unit delay maintaining the polarization. 

14. PMD compensator in optical fiber communication systems applying the method 
in accordance with any one of the above claims and comprising a cascade of adjustable 
optical devices over which passes an optical signal to be compensated and an 
adjustment system which takes the components yl(t) and y2(t) on the two orthogonal 
polarizations from the compensator output signal with the adjustment system 
comprising a controller which on the basis of said components taken computes the 
Stokes parameters S 0 , Si, S2, S3 in a number Q of different frequencies of the 
compensator output signal and emits control signals for at least some of said adjustable 
optical devices so as to make virtually constant the Stokes parameters computed at the 
different frequencies. 

15. Compensator in accordance with claim 14 characterized in that said optical 
devices comprise a polarization controller with control angles fa, fa and two optical 
rotators with rotation angles 0i and 02 and in which said parameters which are adjusted 
consist of said control angles fa, fa and said rotation angles 0i, 0 2 . 



WO 03/050984 



PCT/IB02/05446 



26 



16. Compensator in accordance with claim 15 characterized in that between the 
controller and an optical rotator and between optical rotators there are fibers which 
introduce a predetermined differential unit delay maintaining the polarization. 
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